Velocity gradients statistics along particle trajectories in turbulent flows: 
the refined similarity hypothesis in the Lagrangian frame 
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We present an investigation of the statistics of velocity gradient related quantities, in particluar 
energy dissipation rate and enstrophy, along the trajectories of fluid tracers and of heavy/light par- 
ticles advected by a homogeneous and isotropic turbulent flow. The Refined Similarity Hypothesis 
(RSH) proposed by Kolmogorov and Oboukhov in 1962 is rephrased in the Lagrangian context and 
then tested along the particle trajectories. The study is performed on state-of-the-art numerical 
data resulting from numerical simulations up to Re\ ~ 400 with 2048 3 collocation points. When 
particles have small inertia, we show that the Lagrangian formulation of the RSH is well verified for 
time lags larger than the typical response time r p of the particle. In contrast, in the large inertia 
limit when the particle response time approaches the integral-time-scale of the flow, particles behave 
nearly ballistic, and the Eulerian formulation of RSH holds in the inertial-range. 

PACS numbers: 47.27.-i, 47.10.-g 



I. INTRODUCTION 

One of the most prominent features of turbulent flows 
is the strong variation present in the energy dissipa- 
tion field, a phenomenon called intermittency In an 
attempt to describe quantitatively intermittent fluctu- 
ations in the inertial range of turbulence, Kolmogorov 
and Oboukhov in 1962 proposed a general relation 
linking velocity fluctuations, measured at a given spatial 
increment S r u = u(x + r,t) — u{x,t), with the statisti- 
cal properties of the coarse grained energy dissipation, 
c r = r~ 3 f A , r s e(x, t) d 3 x averaged over a volume, A(r), 
of typical linear size r: 



(1) 



where ~ means "scales as" or "equal in law". Equa- 
tion ([T]) is known as the Refined (Kolmogorov) Similar- 
ity Hypothesis (RSH) and it is considered to be one of 
the most remarkable relations between turbulent velocity 
fluctuations: Many efforts in the last decades have been 
devoted to its validation ([1, 0, @). The importance of 
RSH cannot be underestimated: it bridges inertial-range 
properties with small-scale properties, supporting the ex- 
istence of an energy cascade mechanism, statistically lo- 
cal in space. So far, a rather strong evidence supports 
the validity of the RSH in the Eulerian frame (i.e. the 
laboratory frame). On the other hand, no investigation 
has been reported in the literature on the validity of RSH 
in the Lagrangian frame (i.e. along fluid particle trajec- 
tories). The main difficulty in studying RSH in a moving 



reference frame stems from the necessity to make multi- 
point measurements along particle trajectories in order 
to calculate the stress tensor. As a result, no experi- 
mental measurements along particle trajectories of ve- 
locity gradients exists for time long enough to be able 
to evaluate temporal correlations. Also numerical ex- 
periments are very demanding, requiring refined compu- 
tations of velocity differences along particle trajectories. 
This is usually implemented by computing the velocity 
gradients matrix in Fourier space, then transforming it to 
physical space by (inverse) Fast Fourier Transform, and 
performing off-grid interpolations of the gradients at the 
particle positions. Here we report the first of such mea- 
surements using high-resolution Direct Numerical Sim- 
ulations (DNS) investigations. We also note that when 
the particles transported in a turbulent environment have 
non-negligible size or mass, i.e. they are inertial paricles, 
their trajectories becomes strongly sensitive to the sta- 
tistical and topological properties of the advecting flow 
(0 H i, G3)- The possible validity of Lagrangian RSH 
in this context is far from being trivial and may shed 
new light on the physics of particulate transport in tur- 
bulent flows: an ubiquitous phenomena in nature and in 
industrial applications alike. 

In the present study we will extend the RSH relation to 
the temporal domain and test its validity along the tra- 
jectories of fluid tracers and of inertial particles whose 
density is smaller /larger than the fluid one while their 
sizes span the interval from the dissipative to the inertial 
range of scales. The manuscript is organized as follows: 
First we gives details on the numerical methods of the 
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DNS. Wc then present the extension of RSH to the La- 
grangian domain and we test it on the trajectories of 
fluid tracers. In the last section we investigate the case 
of inertial particles: we show under which conditions the 
Lagrangian RSH still holds and how it should be modified 
in the special case of highly-inertial particles. 



II. NUMERICAL METHODS 

The incompressible fluid velocity u(x,t), V • u = 0, 
evolves according to the Navier-Stokes equations : 



Du du Vp 
— — = — — h u ■ Vu = 

Dt dt p f 



vAu + f. 



(2) 



where p denotes the pressure, pj the fluid density as- 
sumed constant, and / an external large-scale forcing 
injecting energy at a mean rate (e) = (u ■ f). Together 
with the Eulerian field we integrated the Lagrangian evo- 
lution of fluid tracers: dx(i)/dt = v = u(x(t),t), and 
point-particles by means of a model of dilute suspensions 
of small pas sively advected spherical particles, as derived 
in Refs. [UEl El: 



dx 

"dT 



dv „ Du 1 , 
dt Dt t„ 



where x and v denote the particle position and velocity, 
respectively. In Eq. (|3|) the coefficient j3 = 3p//(p/ + 
2p p ) is related to the ratio between the density of the 
particle (p p ) and of the fluid (p/); t p = a 2 /(3(3v) is the 
particle response time, with a the particle-radius. The 
Stokes number of the particle is defined as St = t p /t v , 
where r r) = (v / (e))^ 1 ^ 2 is the dissipative time-scale of the 
turbulent flow. In our simulation the parameters (3 and 
St can be varied independently, therefore it is possible to 
consider also the case (ft = Q,St > 0), corresponding to 
the limit to very heavy particles for which the fluid added 
mass is negligible while Stokes drag is the only relevant 
dynamical force. On the other hand the situation (/? = 
l,St = 0) is equivalent to the case of a perfect fluid 
tracer. 

Eq. @ is numerically integrated by means of a stan- 
dard internally 2/3 de-aliased pseudo-spectral algorithm 
with a second order Adams-Bashforth time-advancing 
scheme. The very same time-scheme is used to track 
the particles evolving according to eq. (J3]) : the time-step 
size in both cases is O(10~ 2 t v ), however particle infor- 
mations are recorded for post-processing/analysis at a 
rate of lO -1 ^. Interpolations of the velocity field, ac- 
celeration field (necessary for ((3J) and velocity gradients 
at the particle positions, are done via a tri-linear algo- 
rithm. For a validation of our numerical method we ad- 
dress Ref. [HI , where a satisfactory comparison on accel- 
eration lagrangian statistics has been performed against 
an independent numerical implementation with several 
different features (field interpolation based on tri-cubic 



scheme, external dealiasing procedure, a slightly differ- 
ent large scale forcing). In our DNS energy is injected 
at large-scale by maintaining the spectral content of the 
first two shells in Fourier space constant. Here we will 
report data coming from two sets of simulations with 
N 3 = 2048 3 and N 3 = 512 3 collocation points, corre- 
sponding to Rca = 400 and 180, respectively, and sam- 
pling the parameter space (3 £ [0 : 3], St € [0:4] with 
64 (f3, St) particles types. A total amount of ~ 10 s par- 
ticles are tracked in time. Results on the clustering of 
these particles in the turbulence have already been re- 
ported in Ref. @, @]. Inertia requires some time before 
particles reach their fractal (or multifractal) statistically 
stationary distributions [H, El- We therefore waited 
till the Lagrangian statistics became stationary (approx- 
imately one large-eddy turnover-time) before performing 
the analysis presented here. Measurements of velocity 
differences and gradients are based on sets of O(10 6 - 10 7 ) 
particles which have been followed in time for few 0(1) 
large-eddy-turnover times. 



III. REFINED SIMILARITY HYPOTHESIS IN 
THE LAGRANGIAN FRAME 



Inertia effect on the statistics of principal 
invariants of velocity gradient tensor 



We have already noted that the effects of inertia may 
be of particular interest for the present study. Inertial 
particles are not distributed homogeneously in the vol- 
ume, centrifugal force tends to concentrate light par- 
ticle inside strong elliptical regions, with high vortic- 
ity 0, H, E3]; and heavy particle in hyperbolic regions, 
typical of intense shear. Following Chong et al. [HI], 
the flow topology may be locally defined in terms of 
the two principal invariants of the velocity gradient ten- 
sor A = Aij = diUj, namely Q = —Tr{A 2 }/2 and 
R = — Tr[A 3 ]/3, (see also [l9| for a recent study). Q rep- 
resents the difference between a rotation-dominated and 
a dissipation-dominated flow topology, e.g., it is positive 
in a vortex core, while negative in a region characterized 
by high strain. The second parameter, R, analogously 
represents the competition between the vorticity produc- 
tion and the dissipation production. Also, the separa- 
trix curve (R/2) 2 + (Q/3) 3 = (so called Vieillefosse 
line [2oT ] ) discriminate between three real or one real and 
two complex-conjugate eigenvalues for A, again meaning 
only-strain or vortical regions. In fig. [T]we show the joint 
probability density function V(Q,R) for different particle 
types as measured in the simulations at Re\ = 180. The 
most striking effect is for light particles ((3 = 3), contrary 
to tracers and heavy particles ((3 = 0) they spend essen- 
tially all of their time in upper half-plane Q > 0, meaning 
that they constantly trapped in vortical regions. 
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FIG. 1: Joint probability density function V(Q*,R*) of Q* = Q/(Q 2 ) 1/2 and R* = R/{Q 2 ) 3/4 , for particles of different types. 
Contour lines are drawn at values 10~ z with z = 0,1, 2, 3, 4, 5, 6 ( from the center to the outside of the figure). The thick line 
traces the curve: (R/2) 2 + (Q/3) 3 = (Viellefosse line), discriminating between complex (above) and real (below) eigenvalues 
of A. Data come from Rex = 180 calculations. 



B. Time correlation of symmetric and 
antisymmetric component of velocity gradient tensor 

One also expects pretty different temporal correla- 
tions between particle trajectories and the underlying 
topology of the carrier flow. We look now at the sym- 
metric/antisymmetric component of the velocity gra- 
dient tensor A, because of their direct link with en- 
ergy dissipation and enstrophy, [2l], [22|, [23|. We show 
in Fig. © the autocorrelation function of enstrophy, 
ft = | (A - A T ) 2 = § V ,,/■ ().,!, - djUi) 2 and energy dis- 
sipation, e = %(A + A T ) 2 = § Y,i.j(9iUj + djUi) 2 , along 
the particle trajectories for different values of inertia. As 
one can see both these quantities have short autocorrela- 
tion time, at Re\ = 180 we find Tq = Cq(t)cIt ~ 7r. n 
and similarly T e ~ 5r r) . However, the autocorrelation of 
enstrophy turns out to be rather sensitive to the type of 
particles, while energy dissipation is probed more or less 
uniformly. This is a clear indication that due to inertia, 
particle tends to leave in regions with very different vor- 
ticity contents, while energy dissipation - although differ- 
ent in intensity - turns out to be a more robust quantity 
in term of coherence-in-timc: This result will be very 
useful in our following discussion. 



C. RSH and its generalized formulation 

Along the trajectory of a fluid tracer x(t) the velocity 
difference will be denoted as 5 T v = v(t + r) — v(t) and 
similarly we define a coarse grained energy dissipation 
measured along the trajectory as e T = f* +T e(x(t), t) dt 



(see also [24j ) . The RSH (flj can be translated from space 
to time by making the assumptions that S T v ~ S r u and 
e T ~ e r when r and r are linked trough the eddy turnover 




FIG. 2: Temporal autocorrelation function, i.e., Cxij) = 
(X'(t) - X'{t + r))/{X' 2 ) with X'(t) = X(t) - (X), of 
the enstrophy X — fl (top) and the energy-dissipation-rate 
X — e (bottom) for fluid tracers and for inertial particles 
with St = 0.5, j3 = 0, 3, at Re x = 180. 



time definition, r(r) ~ r/8 r u. This argument leads to the 
Lagrangian refined similarity hypothesis (LRSH): 



(4) 



In order to test Eq. Q one should verify, for any expo- 
nent p, the scaling relations: 



(5) 



where ~ means equal apart from a multiplicative con- 
stant depending only on p, in the inertial range. In the 
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time domain the inertial range is denned as the inter- 
val, t v -C t <C Tl, where Tl is the Lagrangian inte- 
gral time scale, which is estimated as the autocorrelation 
time of velocity of fluid tracers, i.e., Tl = J °° C v (t) dr, 
with C v (t) = (v(t)v(t + t))/(v 2 ). As one can estimate 
Tl/tti ~ Rex, the extension of the inertial range in dissi- 
pative time-scale units extends over roughly two decades 
in the present numerical study. 

In contrast to the 4/5-law (consequence of the 
Karman-Howarth equation) leading to exact scaling 
properties for third order velocity increments in the Eu- 
lerian frame, we do not have any exact scaling relation 
derivable from NS equations in the Lagrangian domain. 
Furthermore, it is known that in the Lagrangian frame, 
finite Reynolds effects induce larger deviations from a 
power law regime than what observed in Eulerian frame 
PH - To overcome these effects, following (26|, we can 
generalize the above expression ([5]) by using its Extended 
Self Similarity (ESS) form, namely: 

((M p )-( ii gp) P/2 (4 /2 > (6) 



D. Numerical tests of LRSH 

In Figure [3]Ja) we present a test of Eq. ([S} for p = 4 
for particles with 8 = 1, St = 0, i.e. fluid tracers (cir- 
cles) and very-heavy particles with = 0, St = 2 (trian- 
gles). In Figure [3jb), we show instead the relation from 
Eq. ^ for the same particle types. Two major results 
emerge. First the LRSH, as expressed by Eq. (|6|) is well 
verified for the transport of particles in turbulent flows. 
The use of the ESS version for LRSH is able to over- 
come finite size/time effects which are usually observed 
at relatively low Reynolds number (see [H|]). The second 
important result, which will be investigated later on in 
this manuscript, comes from inspecting the validity of ([6]) 
for different Stokes numbers. Equation ([5]) is supposed 
to be valid both in the inertial range and in the dissi- 
pative range (where the velocity field is smooth) though 
with different offset. This is clearly observed in Figure 
|3{b) for the case St = 0. It is already known that by 
increasing the Stokes number, particles tends to escape 
from strong vorticity region, thus decreasing the effect of 
the dip present in between dissipative and inertial scales 
p7| . As a consequence, for St = 2 we observe almost no 
deviation of Eq. ^ in the range of scales between the 
inertial and the dissipative ranges. 

To have a more quantitative check, we look now at the 
ratio between the two sides of Eq. ©, namely at ((5 T v) p ) 
divided by (e p T /2 )((S T v) 2 )P/ 2 (e T )-P/ 2 , as a function of the 
time difference r. In FigureQ]we show its behavior for the 
order p = 4 fluid tracers particles, the time difference r is 
normalized by the dissipative time-scale r^. We observe 
a plateau (see circles symbols, in Fig[3} for t/t v > 5. No- 
tice that also in the dissipative range the compensation 
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FIG. 3: Test of LRSH along the trajectories of tracers and 
heavy particles at Re^ = 400. (a) For p — 4 we show Eq. © 
for St = (circles) and St = 2 (squares), (b) We show the 
validity of Eq. (j6j) for the same values of p and St. Straight 
lines correspond to the theoretical scaling prediction. Data 
come from Re\ = 400 calculations. 




FIG. 4: Test of LRSH along the trajectories of tracers (Re^ = 
400). It is plotted ((S T v) 4 )/( ((<M 2 ) 2 (X 2 T ){X T )- 2 ): (cir- 
cles) represent the case X = e, (squares) the case X = Q, 
and (triangles) the case X = const. In the inset the joint pdf: 
P{e' T ,\5 T v\') and P(Q.' T ,\S T v\') at r = I3r^; (note that the 
prime symbol denotes variables normalized respect to their 
mean values, i.e. x' = x/{x)). 



works well, as it should from the requirement that the ve- 
locity field becomes differentiable, 5 T v ~ r. However, the 
plotted ratio shows a mismatch with the value attained 
in the inertial range. The transition between the two 
plateaux occurs around the dissipative time scale, where 
the presence of vortex trapping has been shown to spoil 
the scaling behavior of Lagrangian structure functions 
{(S T v) p ) of the tracers [H, [23, H3, H3 ■ In the same figure 
we show that using the coarse grained enstrophy, i.e., fl T 
instead of e T , the compensation is worse (squares). Simi- 
larly, compensation with enstrophy does not work neither 
for heavy nor for light particles (not shown). Compen- 
sating without coarse grained quantities, i.e. checking 
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FIG. 5: Same as in Fig|4] here also along the tra- 
jectories of heavy particles (Re a = 400). It is plot- 
ted {{8 T vf){e T ) 2 / ({e 2 T ){{8 T vf) 2 ). Particles with St = 
0.6,1,2,10 are compared with the result for tracers (solid 
black line). The LRSH is satisfied both in the inertial and 
in the dissipative ranges. The prefactors, Ai and Ad, how- 
ever differs in the two regions. Notice that for the largest 
Stokes, St = 10, the smallest time lags where LRSH is veri- 
fied, as expected, roughly r ~ lOr,,. In the inset, the behavior 
of ratio of prefactors Ad/Ai is plotted vs. St. For small St 
values a behavior as St~ ' 38 is found (solid line in the inset). 

the deviation from dimensional, non- intermittent, scal- 
ing does not provide a good plateau (triangles in FigSJ). 
This result supports the validity of LRSH only when us- 
ing the energy dissipation as the main driving process 
along the particle motion. The behavior for intense fluc- 
tuations (moments higher then p — 6) can not be checked 
quantitatively due to the lack of statistics. Nevertheless, 
in the same figure, we show the joint probability density 
functions, V(fl T , \$tv\) and V(e T , \d T v\), for a time lag 
t = 13 Tq. Velocity increments are more correlated with 
coarse-grained energy than with enstrophy, as shown by 
the high probability measured for simultaneous intense 
values of \S T v\ and e T . 

Having established the validity of the LRSH, we strive 
now at investigating the effect of different Stokes number 
and different density properties. 

IV. LRSH IN THE (J3, St) PARTICLE 
PARAMETER SPACE 

A. Heavy particles (0 = 0) at St ~ O(l) 

When particles have inertia their trajectories deviate 
from material lines of the flow. In principle, one expects 
that for very small value of the inertia (particles very 
close to fluid tracers) no appreciable discrepancies can be 
measured. In Figure [5] it is shown the test for the LRSH 
compensated with the energy dissipation rate. From One 
can appreciate that in the inertial range, e.g. t/t^ > 5 
the LRSH is well verified for all the Stokes considered. In 



the dissipative range it is also verified but with a different 
proportionality constant. In particular, the important 
mismatch observed between the two plateaux for tracers 
in (FigH]) here is reduced considerable as soon as some 
inertia is switched on. This confirms that heavy particles 
arc quickly expelled out of vortex filaments, and there- 
fore much less sensitive to the transition around r / t v ~ 1 
than tracers [13] (the opposite will happen for light par- 
ticles, see below). 

The behavior of the ratio (Ad/Ai) of the plateaus dis- 
played by ((<5 T u) 4 )(e r ) 2 / ((s 2 ) ((S T v) 2 ) 2 ) respectively in 
the dissipative range (Ad for r <C t v ) and in the inertial 
range (Ai for r < 10t v ), is shown in the inset of Fig. [5] 
The estimate for the slope of Ad /Ai vs. St can be pro- 
vided by the following reasoning. First we notice that the 
inertial constant Ai is almost insensitive from the Stokes 
value, therefore the dissipative constant Ad carries all 
the St dependency. Moreover, we have measured that the 
single point energy dissipation statistics is pretty insensi- 
tive to the Stokes number (see again fig. ^ . As a conse- 
quence the main dependency on St for the ratio Ad/Ai 
comes from the flatness factor F(t) = ({Stv) 4 ) / ({5 T v) 2 } 2 
in the intermediate-dissipative t limit. It is reasonable to 
estimate the difference between F(t) at changing Stokes, 
but fixed Reynolds, as given by the value of the flatness 
at the particle response time: F(t p ) ~ St l > 4 -~ 2 '' 2 , where £ p 
is the p-th. order scaling exponent for Lagrangian struc- 
ture functions ((5 T i>) 4 ) *~ r* p . Based on the experimen- 
tal values ^4 — 1.6 and (2 ~ 1 [30| . this estimate gives 
A d /Ai ~ St-° A , not too far from the fit Sft-o.38io.05 to 
our numerical data, see inset of Fig [5j 



B. Finite density contrast: heavy and light 

particles 

We now look at the statistical properties of particles 
with finite density contrast, i.e. also j3 ^ 0. In Figure [S] 
it is shown, for St = 1, the behavior of the compensated 
tests for LRSH for different values of j3 spanning the full 
range [0 : 3]. In the inset it is also shown the behav- 
ior, this time as function of (3, of the Ad /A] ratio. No- 
tice how the critical value (3=1, discriminating between 
heavy (/3 < 1) and light (/3 > 1) particles plays a crucial 
role. Again LRSH is well verified in the inertial range, 
but the change to a different plateau around t/t v ~ 1 is 
now much more abrupt when light particles are consid- 
ered: for those the vortex trapping is more pronounced, 
as all light particles quickly move towards high vorticity 
regions, showing a very sensitive dependency around the 
dissipative time dynamics. A model for the dependency 
of Ad/Ai vs. (3 is presently not available. 

C. Heavy particles (/? = 0) with large inertia 

Having studied the case of particles with small inertia, 
we now focus on the case of extreme inertia, i.e., when 
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FIG. 6: Same as in fig. [4]along the trajectories of heavy/light 
inertial particles with St — 1.0 (Re a = 180). It is 
plotted ({S T v) 4 )(e T } 2 / ((el){(S T v) 2 ) 2 ). Particles with j3 = 
0,0.5, 0.75, 1.0, 1.5, 3.0 are compared with the result for trac- 
ers (solid black line). The LRSH is satisfied both in the in- 
ertial and in the dissipative ranges. As for the case of heavy 
particles the prefactors differs in the two regions. In the inset, 
the behavior of ratio of prefactors Ad/Ai is plotted vs. 0. 




FIG. 7: Particles with very large inertia do not verify the 
LRSH (HU, but follow the Eulerian version (0. We show this 
for the order p — 6 on particle trajectories with St = 70, 
at Rex = 400. As it can be seen, in the inertial range, the 
Eulerian RSH compensate better than LRSH. 



the response time of the particle is at the top end of the 
inertial time-range, or St ~ O(10). In this condition, 
for time lags r < r p when the particle filters out most 
of the underlying turbulent fluid fluctuations and evolves 
nearly ballistically, one can predict a different behavior 
for 5 T v. Along the trajectory of a ballistic particle, the 
relation linking scale to time is r(r) ~ r/vo where the 
typical particle velocity vq is proportional to root mean 
square fluid velocity. Recasting eq. JT]) from space to time 
notation we obtain again an Eulcrian-likc RSH relation, 



fx \ II \V3 1/3 

(5 T v) ~ {t/vo) £ T , or: 

((6^}* (i.y 3 {E*/ 3 ) (?) 

The generalized version of ([7]) reads now 

/-J>/ 3 \ 

((^-^-^((M) 3 }^ 3 (8) 

In Figure [7| we present a test of this idea for ((6 T v) p ) 
with p = 6 : For particles with very large Stokes numbers 
(St = 70) we compensate the velocity increments both 
with respect to the prediction of the Lagrangian RSH and 
with respect to the prediction of the Eulerian RSH in its 
generalized version. The compensation with the Eulerian 
RSH works appreciably better in the range r < St ■ t v 
than the compensation with LRSH. 

V. CONCLUSIONS 

In summary, some important statistical properties 
of velocity gradients along trajectories of fluid tracers, 
heavy and light particles have been investigated. We used 
high-resolution, high-statistic numerical data to correlate 
the temporal properties of velocity gradients and veloc- 
ity differences along trajectories. We demonstrated that 
the Refined Similarity Hypothesis is well verified both 
for fluid particles and particles with response time in the 
dissipative regime, a feature that we dubbed Lagrangian 
RSH. Around the dissipative time lags, heavy and light 
particles behave strongly differently, due to the effect 
of being expelled/concentrated out/in vortex filaments. 
The dynamics at those time lags becomes markedly de- 
pendent on the underlying topological flow properties. 

Understanding the RSH in the Lagrangian domain may 
also have important applied consequences. In many ap- 
plications, the geometry of the system and/or the in- 
tensity of turbulence do not allow for a direct attack of 
the problem using numerical simulations of the Navier- 
Stokcs equations. Modeling is needed for both the un- 
derlying fluid and for the particle equations. Typically, 
the ideal model, would like to replace Eqs.Q-© with a 
Langevin-like equation for the particle evolution [32l . |33| : 
dx/dt = v, dv/dt = V(A)v + T(t) where the V rep- 
resents some stochastic noise induced by the underlying 
turbulent fluctuations. The hard physical problem is in 
the modelization of the drift term, 2? (.4), depending on 
the loc al g radient structure along the trajectories (see 
[H, [35|, [36[ for recent attempts) . Such term should also 
take into account effects induced by preferential concen- 
tration in/out vortex filaments for the case of inertial 
particles around the dissipative time lags. 
The numerical database presented here can play a crucial 
role for benchmarking stochastic models for tracers and 
inertial particles in turbulence. Data from this study are 
publicly available in unprocessed raw format from the 
iCFDdatabase ( |http: //cf d. cineca. it[ ). 
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During the preparation of this manuscript we got aware 
of a slightly similar investigation [37| where Lagrangian 
correlation of velocity and pressure gradients are studied 
conditioning on the initial Eulerian energy dissipation, 
a sort of mixed Eulerian-Lagrangian refined Kolmogorov 
hypothesis, different from the fully Lagrangian view point 
adopted here. 
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